MPI-PhT/2001-14 
hep-ph/OTU^ 

Neutrino anomalies and large extra dimensions 



Amol S. Dighe^* and Anjan S. Joshipura^^ 

^ Max- Planck- Institut fur Physik, Fdhringer Ring 6, D-80805 Miinchen, Germany. 
^ Theoretical Physics Group, Physical Research Laboratory, 
Navarangpura, Ahmedahad, 380 009, India. 



Abstract 



Theories with large extra dimensions can generate small neutrino masses 
when the standard model neutrinos are coupled to singlet fermions propagat- 
ing in higher dimensions. The couplings can also generate mass splittings and 
mixings among the flavour neutrinos in the brane. We systematically study 
the minimal scenario involving only one singlet bulk fermion coupling weakly 
to the flavour neutrinos. We explore the neutrino mass structures in the brane 
that can potentially account for the atmospheric, solar and LSND anomalies 
simultaneously in a natural way. We demonstrate that in the absence of a 
priori mixings among the SM neutrinos, it is not possible to reconcile all these 
anomalies. The presence of some structure in the mass matrix of the SM neu- 
trinos can solve this problem. This is exemplified by the Zee model, which 
when embedded in extra dimensions in a minimal way can account for all the 
neutrino anomalies. 
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I. INTRODUCTION 



The present data from the experiments on atmospheric, solar and reactor neutrinos 
indicate neutrino flavour oscillations. The data from each of these sets of experiments 
individually can be explained by a single dominating mass square difference Am? and a 
mixing angle 9. The atmospheric neutrino data |l| indicate ^ as the dominant mode, 
with Am^j^ = (1 — 8) ■ 10~^ eV^ , sin^ 2Q = 0.8 — 1.0. There is no compelling evidence that 
the electron neutrinos participate in the oscillations of atmospheric neutrinos. Moreover, the 
CHOOZ experiment 0] gives an upper bound on the mixing of z/g: we have sin^ 26'e < 0.1 
where z/3 is the mass eigenstate such that Amg^^ ^ Am^g ~ Am^^^. The three MSW 
solutions (LMA, SMA and LOW) as well as the vacuum oscillations can provide reasonable 
fits to the solar neutrino data 0], all these solutions have AttIq < 2- 10~^ eV^. The results of 
the LSND experiment are neither confirmed nor fully excluded by the KARMEN2 data [Q, 
and the combined fit allows a region Am^^.^^, = (0.1 — 1) eV^ , sin^ 29^e ~ 10~^ — 10^^. 

In the context of only three known neutrino species, the Am^s corresponding to the so- 
lutions of the three neutrino anomalies (atmospheric, solar and LSND) cannot be reconciled 
and if the oscillation mechanism is to be used to explain all the anomalies^, the introduction 
of a sterile neutrino becomes necessary. There are two simple schemes of mixing among four 
neutrinos which can account for all the neutrino anomalies: 

(i) In the "3+1" scheme |^, the masses of the three active states are separated from that 
of the sterile one by the LSND scale. The possible — oscillations seen at LSND arise 
through simultaneous mixings of Ve and with the sterile neutrino z/^. Such a picture is 
strongly constrained by laboratory disappearance experiments. It was argued ^ that these 
experiments constrain the Ve — mixing below the experimental signal at LSND. This has 
changed with the new LSND data and the 3+1 scheme is claimed P| to be viable for some 
parameter range, though the detailed statistical analysis in [|10| does not support this claim. 

(ii) The second scheme — known as the "2+2" scheme [§|ll| — corresponds to two nearly 
degenerate pairs of neutrinos separated by the LSND mass scale. The mixings within the 
pairs account for the solar and atmospheric flux deficits. The small oscillation probability 
seen in the LSND experiment can be explained by requiring that z/g and belong to two 
different degenerate pairs. This implies that the sterile state has significant admixture 
either with u^, or z/^ and thus plays an important role in generating either the solar or the 
atmospheric neutrino deficit. The data on atmospheric neutrinos disfavour z/^ — z/^ as the 
dominant source of the z/^ oscillations |jl|. The global fit to the available sets of the solar 
neutrino data are also better described in terms of Vp conversion to active rather than the 



sterile state, although the latter possibility is not completely ruled out [|T2l. The 2+2 scheme 
will be ruled out if both the solar and the atmospheric neutrinos are found to be mainly 
oscillating to active neutrinos. 

It follows that both the above schemes are strongly constrained and may not prove 
adequate for incorporating the LSND signal. In this case, the description of data in terms 
of four neutrino mixing would be challenged and more complicated possibilities would be 



^Some non-oscillation "exotic" solutions have been proposed ||^, but we shall not discuss them 
here. 
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needed. Quite apart from phenomenological incompatibility, the very existence of a light 
sterile state needs theoretical justification. While there have been attempts in this direction 
I3| , the introduction of more than one light sterile states can add to theoretical difficulties. 
The above theoretical and phenomenological problems find natural solutions in theories 



with extra dimensions of mm size |T^. Firstly, such theories can provide justification for the 
existence of light sterile states which form a part of the Kaluza Klein (KK) tower associated 
with a singlet fermion residing in the bulk volume of the full higher dimensional theory. In 
addition, the presence of more than one light sterile state allows additional phenomenological 
possibilities not envisaged in the context of the 2+2 or 3+1 schemes mentioned above. As a 
result, these theories provide a potentially useful framework for understanding all neutrino 
anomalies. 

The standard model particles are assumed to be confined to a three-brane embedded in a 
higher dimensional bulk. The gauge singlet states are not confined this way and propagate in 
the bulk. Such theories have been argued to provide some understanding of neutrino masses 
if one assumes the presence of a singlet neutrino ub propagating in the bulk. The couplings 
of ub to the fiavour states residing in bulk are suppressed by a volume factor corresponding 



to the extra dimensions and lead to very small neutrino masses p^JT6[| in spite of the new 
physics at TeV scale. 

The minimal approach consists of assuming only five dimensions^ and the presence of 
only one bulk neutrino coupling to all the known fiavors through 

— la(pl^B + h.C. , (1) 
V 

where 



rrir, 



Mp V TeV 



6-10-^eV/iaUTT7 • (2) 



Here ha is the Yukawa coupling, and Mf[Mp) is the fundamental (Planck) mass scale. 
Neutrino masses and mixings are determined by four parameters: three mixing parameters 
= V2maR and an overall mass scale = ml + + m^. 

From the point of view of the 4-dimensional theory, the ub residing in a five dimensional 
world contains a massless state which can provide a sterile neutrino needed to understand 
various neutrino anomalies. The n^^ KK excitation associated with ub has a mass m„ = n/R. 
For R ~ 0.1 mm, masses of some of the low lying excitations would be in the range needed 
for the occurrence of the MSW resonance for solar neutrinos. The mixing of the KK states 
with z/e is governed by ^e- For Mf in the TeV range, this mixing is in the range needed for 
the small mixing angle solution of the solar neutrino problem. Hence, the MSW solution in 



this context is correlated to observable extra dimensions and TeV mass scale |17]. However, 



these values of the parameters cannot provide the atmospheric neutrino scale if Yukawa 



^ The observed accuracy of Newton's law over earth-sun distance requires the presence of two or 
more extra dimensions. We have implicitly assumed that the radii of all other extra dimensions 
are much smaller than that of the fifth dimension. One can work in this case with an effective five 
dimensional theory. 
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couplings ha < 1- This makes it difficult to reconcile all the neutrino anomalies and one 
needs to go beyond this minimal picture. The following different possibilities have been 
considered in the literature: 

(i) Extension of the minimal scheme with three bulk neutrinos was considered in 



The phenomenology of this case is described in terms of an arbitrary mixing matrix, three 
Dirac brane-bulk masses rria and corresponding Ua — mixing parameters E,a- The solar 
and atmospheric neutrino deficit can be explained for small S,a- But this requires larger 
than the natural value expected from (|]), and consequently the fundamental scale of 
about three to four orders of magnitude higher than TeV. On the other hand, for larger ^2 
and/or ,^3, the atmospheric neutrinos oscillate to the unfavoured mode of the sterile state. 
Moreover, one cannot explain the LSND results in a satisfactory manner in any of these 
cases W^- 



(ii) Lukas et al. [|19| considered an extension which allows for a mass term of the bulk 
neutrino. Such an extension is compatible with all neutrino anomalies provided one intro- 
duces three bulk neutrinos and lepton number violation in the brane-bulk couplings. Some 
of these couplings are required to be as large as ~ 10~^ eV as in the case (z). Thus either 
one needs a relatively large fundamental scale or Yukawa couplings much larger than 1. 

[iii) The neutrino masses and mixings exclusively occur due to physics in the bulk in 
both the above cases. It is conceivable that physics in the brane itself may have some seed 
for the neutrino oscillations to occur. This possibility was considered in It was shown 



that if the brane neutrinos have some masses then their coupling to the bulk neutrino can 
induce mixings, and hence oscillations, among the brane neutrinos. In this context, models 



have been considered |]2T| that try to solve some of the neutrino anomalies. 

The aim of this paper is to study the feasibility of the suggestion {iii) from the point of 
view of solving all neutrino anomalies. Motivated by the observation that strong mixing of 
sterile state with is disfavoured [jll, we shall work in the limit of only small brane-bulk 
mixing, i.e. small E,a- Neutrino oscillations can be largely described in this framework from 
an effective 4x4 matrix in the space of three active neutrinos and the massless mode of 
vb- One can use the standard seesaw approximation to determine this matrix provided one 
allows non-leading terms in this approximation. We present the general formalism for this 
in Sec. ||. Non-zero brane mass for the electron neutrino can significantly influence the 
solution to the solar neutrino problem proposed in |]T^. We discuss this effect quantitatively 



in Sec. |T|. These results are then used in Sec. |V| to study phenomenology of the minimal 
case in which one allows neutrino masses but no mixing in the brane. It is shown that 
one cannot solve all the neutrino anomalies in this case. In Sec. 0, we study a specific 
structure of the neutrino mass matrix based on the Zee model |^ and show that this, when 
embedded in the extra dimensions, can explain all the neutrino anomalies. The last section 
gives summary of our results. 

II. GENERAL FORMALISM 
A. The Lagrangian 

We start by taking only one singlet neutrino in the bulk, with the free bulk action given 

by 
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5buik = / d^x dy tdA^ , (3) 



where 7"^ {A = 0, 1, 2, 3, 4) are the 5-dimensional Gamma matrices. Here x = x^, x^, x^) 
are the usual 4-d coordinates, and the SM brane is at 7/ = = 0. Let the bulk neutrino 
couple to all the known active flavours, so that the brane-bulk Yukawa coupling i^ 

^brane = " / d'x ( L^H^ + h.c] . (4) 

Jy=o \\'2'kRMp J 
The four dimensional effective action may be obtained by expanding ^ as 

'^{x.y) = ^=^T.'^n{x)exp('-^) ' (5) 

where each of the modes may be written in terms of two 4-d left-handed Weyl spinors 
^,77 as 

The effective 4-d Lagrangian involving KK modes can be written as 

- Cmass = ^a^aiYl + ^n^nVn + h.C. (7) 

a neZ neZ 

where = n/R. Defining new fields 

= ± ^-n)/V2 , Vn = iVn ± V-n)/V2 , (8) 

we can write the above Lagrangian (|^ as 

- ^mass = J2 ^aJ^aiVO + ^ J2 Vn) + J2 ^niCVn + CVn) + h-C. (9) 

The fields and rj~ decouple completely from the rest. 

We add to the above Lagrangian a general mass term for the flavour neutrinos (« = 
e,/^,r)g: 



^ Since and ^ have the same transformation properties under 5 dimensional Lorentz symmetry, 
we have chosen to name as ^ the bulk spinors with lepton number +1 (as in [^), so that the 
brane-bulk coupling (Q) is lepton number conserving. Lepton number violating terms of the form 
h'^LaH^^ are possible in principle, but we can get rid of them through the imposition of a Z2 
symmetry [^. The same symmetry also forbids a possible Dirac mass term in (|^. 

This term explicitly violates the lepton number in the brane, but it is still allowed under the Z2 
symmetry that forbids the bulk-brane lepton violating coupling [p^]. 
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The net Lagrangian relevant for the neutrino oscillations can then be written as 

1 



Cmass = -{l^e, Z^/., ^^r, Vo, V, -^v (^^e, ^^/x, ^^r, ?7o, C) 



where 



The Majorana mass matrix A^,^ is 













\ 




IJ'fj.fj. 






V2M^ 


















me 




rrij. 























V 













; 



where the matrices Mq, and M^^ are defined as 
M„ = m„(l,l,l,...) 



iDiag(l,2,3,...) . 



B. Diagonalization of the neutrino mass matrix 

We now use the standard seesaw approximation to diagonalize eg . ([T3|) . This 
written in the block form 



rriL m 
M 



where 



rriL 









me ^ 




( V2Me 


0\ 








m^ 




V2M^ 









f^TT 


rrir 


, m = 


V2Mr 





V me 


mf. 




J 




[ 





The matrix M is diagonalized by the unitary matrix Y: 
Y = ^ ~^ \ so that M = Y^MY 



M = 



M 



Mrii 



-M 



Then the action of the matrix Y = \ ^ I on A^,^ will be 

or 



M 



where m = mY 



( M, -Me\ 
\Mr -Mr) 



(18) 



The matrix M.^ may be "block-diagonalized" in general by 
Vl - BB"^ B 



W 



-fit VI - B^B 



, such that My = W MyW 



fflL 

M 



(19) 



Here, in the seesaw approximation that is valid for fiai3,^a ^ ^/R, we may expand the 
relevant quantities as 



itll = mi — mM m — -{rhM rh mi + h.c.) + ... 

~ ^ 1 

M = M + -{m^mM~^ + h.c.) + ... 
B = rhM~^ + rriLrhM''^ + ... 



(20) 

(21) 
(22) 



In the above expressions, we have retained the non-leading terms neglected in the usual 
seesaw approximation. This is necessary because of the fact that the leading term 



rhM m 







(23) 



vanishes identically as can be checked from eqs. ( [T^ ) and (|TSp. The non-trivial correction 
to rriL is thus provided by the third term in eq.(^) and the is given by 

rfiL = rriL — -{mM~'^m^mL + h.c.) 



rriL - 



I Xl XeXt^ XeXr 0\ 

XeXl-t XfM Xi-iXt 

XeXr X/iXt Xt ^ 

V o; 



rriL + h.c. 



(24) 



where we have defined Xa = \Jj2-^^a = f^aRyYl ^ for future convenience. 

Let K be the matrix that diagonalizes niL through K^ifiiK = Diag(mi, m2, m^) and V 
be the matrix that makes the elements of M positive: 



V 



1 

i 



so that V^MV 



M^g 



M 



(25) 



The matrix A^j, is then completely diagonalized by 
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K 



K 
V 



through K'^MyK = Md 



(26) 



The net mixing matrix which diagonahzes A^j, is then 



U = YWK = 

To the first order in niaR, we have 

/ nif, \ 



Vl - BB^ K 



BV 



-YB^k fVl - BB^ V 



(27) 



B ^ mM~i = R 



rrir 

V ; 



( Mi/„ -Mi/„ ) where Mi/„ = (l 



1 1 

2'3' 



(28) 



Since B ~ mai? ^ 1, we may expand in powers of BB'' to write the flavor states ( 
and = rjo) as 

= ^[5a/3 - ^{BB^)a,p]ki3iUi + ^ B^^niVn + ^C) 
/3 ^ neN 

= E['^-/3 - l{BB^)^f,]ki,,iy, + V2Y^ B^nNn 

j3 ^ neN 



(29) 



where iV„ = (77+ + zC)/V2- Note that 

V2Ban = 



^a/n ae {e,/i,r} 
a = s 



(30) 



so that the mixing between the KK excitations and the z/„, t^q is of 0{^) as in the standard 
seesaw approximation, and 



(BB- 



XaX/3 a, /? e {e, /i, r} 
a = s OT [3 = s 



To the first order in rriaR, the approximation 



U K 



(31) 



(32) 



is vahd for a G {e, /x, r, s} and i G {1, 2, 3, 4}. 

Eq. (^9]) coincides with the corresponding equation in case of a single flavour derived 



for example in ||T^. The mixing of higher KK states with the brane states is unaffected 
by the presence of masses in the brane. In contrast, the mixing among four neutrinos as 
determined by Kpi in (^) is strongly dependent on the presence of neutrino masses in the 
brane. In particular, the mixing among brane neutrinos may be generated entirely from 
their couplings to the bulk neutrino. We shall study implications of this in Sec 
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III. ELECTRON NEUTRINO MASS AND SOLAR ANOMALY 



The conversion of Ve from the sun to sterile neutrinos can occur in this theory through 
the resonance of v^. with the tower of the KK states with masses ^. This has been explored 
in detail [jl^ in the case of a massless electron neutrino. With R ~ 10^ eV~^, the KK 



states are separated by mass differences of the order of 10~^ eV, so the v^, mass A < 10"^ eV 
will still not affect the v^. survival probability calculated therein. However, in some of the 
scenarios that we will be considering in this paper, z/g (more precisely, the mass eigenstate 
with a considerable fraction of the electron flavour) can have a non-zero mass A > 10~^eV. 



We discuss in this section modifications introduced in the treatment of [T^ due to such a 
mass. In particular, relatively large v^. mass can spoil the solar neutrino solution proposed 
therein. 

For an electron neutrino with energy the density of the layer of resonance with the 
n**" KK state is 

^""U^ ^ j 2EG^(Fe - r„/2) ' ^^^^ 

where mjq is the nucleon mass, Gp is the Fermi coupling constant, and Ye{n) is the number 
of electrons (neutrons) per nucleon in the medium. Eq.(^) gives the set of energies EnR 
beyond which Vf, undergoes resonance with the n*'* KK state: 

E^R ^ 0.3 MeV ■ {n^ - \^B?) ( j , (34) 

The values of EnR can be appropriate for solar neutrino provided R ~ 10^ eV~^. The specific 
excitations which participate in the resonance are influenced by the value of \R. 

The neutrino conversions take place mainly in the resonance layers (PB|). Since the mixing 
angle is small (^e/''^ 1), the resonance layers are well separated. As the neutrinos travel 
outwards from their production point inside the sun, they encounter the resonance layers 
corresponding to the densities ([33|). The survival probability after passing through the 



resonance layer is 



Pn ~ Exp 



A^ 



2\R^ J n^E dp/dr 
The net survival probability is 



p=Pn 



(35) 



surv 

n=n 



n Pn (36) 



where rimax corresponds to the highest density resonance that the neutrinos with energy 
E encounter, and rimin is determined by the condition n > XR. Since P„ approaches 1 
rapidly (-P„ ~ Exp[— C/n^]) with increasing n, the upper cutoff in n is not of much practical 
significance. The lower cutoff however may have a significant impact on the final survival 
probability, depending on the mass. 

For A = 0, the value of P„ is independent of the specific excitation n which takes part 



in the resonance |17|. This gets altered when A is of the order of the MSW scale or larger. 
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In particular, a substantially large value for A implies large n^m for the resonance and the 
suppression of Psurv due to P„s with low values of n is absent. This implies a higher z/g 
survival probability. 

We show the value of Psurv as a function of z/g energy in Fig. |l] for different values of A 
using the parameters 

4^2 _ iQ-3 ^ ji=iq3 gv-i , Ye = Yn = 0.5 , p/{dp/dr) = 1.4 ■ 10^^ eV"^ . 

The value of Psurv is seen to be significantly affected by a non-zero A. In particular, it 
becomes ~ 1 for A ~ eV. In this case, the resonance with the KK states is ineffective 
in converting the solar neutrino with small mixing C,e- We shall encounter different cases 
considered here in the subsequent phenomenological analysis. 



IV. UNMIXED ACTIVE SECTOR 

In this section, we discuss specific forms oirriL which correspond to massive but unmixed 
neutrinos in the brane. Mixing among them can arise only indirectly through their couplings 
to the bulk neutrinos. This possibility was proposed in ||2^ where the neutrino oscillation 



patterns were studied numerically in this scenario. We give below an analytic discussion and 
concentrate on the feasibility of solving all neutrino anomalies in this context. Many of the 
important features needed for this can be elucidated by considering only two generations. 
We thus first consider the case of two active flavours a and (3. 



A. Two generations with Majorana mass 

We take the following form for the flavour mass matrix /Zq,^: 

Pi ^ 

P2 



(37) 



The 3x3 Majorana mass matrix in the basis of the three "flavor states" i/q, z//3, Vs is given 
by 



rriL 



\ 



I Pi rria 
P2 rri/s 



(38) 



After taking into account the mixing with the bulk modes, the mass matrix rfiL becomes 
(see eq.(ll)) 



rriL 



( /Wi(l-Xa) -(/Ui +/i2)XaX/3/2 ^ 

{Pl + P2)XaXp/'^ P2il - Xp) 

mpt / 



niat 



(39) 



where t = 1 — xi/2 with xi. 



Xl + X% 
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As expected, the bulk states have generated mixing among the flavour states. The 
corresponding mixing angle is seen to be given in the limit /ii.2 ^ ^ti^,/? by 

tan 26 = r r 2 ?T ■ (40) 

(/i2 - /il) - (/i2X^ - f^lXl) 

This induced mixing among flavour state can be large if fii ~ /i2 while it tends to be small 
for /il ^ — yU2- This is to be contrasted with the pseudo-Dirac case where large mixing is 
linked to the presence of almost equal and opposite masses. The mixing angle 6af3 is maximal 
in the limit {fii = fi2,Xa = X/s}- The brane neutrinos are degenerate in this case but this 
degeneracy is lifted by couplings to the bulk neutrino. Thus both the mass splitting and the 
(large) mixing can be completely attributed to the higher dimensional physics in this case. 
The mass squared difference among the three neutrinos can be worked out from (^) and 
are given by 

Aml^ = 2fi\-xl + ^tan^2<f)) , Am^, ~ A32 ~ /x^ , (41) 



where 



/i = /il = /i2 , tan2(p=-- ^ , r_i = ^ . (42) 

1 -XI 



If the flavour states {a, (3} are identifled with {i^^, z/^}, the the large mixing in (|40| ) can 
be identifled with the mixing observed in the atmospheric anomaly. The relevant (mass)^ 
difference is given by Ami2 in (^) . A potentially interesting range of parameters is /i ~ 1 eV 
and ~ 10~^. This range reproduces the atmospheric mass scale, has the required large 
mixing, and also has the seed of explaining LSND result through /i ~ eV when the coupling 
to a nearly massless is turned on. The extension of this scheme is discussed in Sees 

Let us consider the alternative possibility corresponding to {a = e, P = fi}. Depending 
on the ranges of parameters, particularly the mass /i, solution of the solar neutrino problem 
in this context has several interesting aspects: 

{i) If fi is chosen to be around the MSW scale or lower (/i^ < 10~^ eV^), then can 
get converted to the KK excitations through MSW resonance with them, as discussed in 
Sec. |T1|. Even if the conversion probability is somewhat reduced when /x ~ MSW scale, such 
conversions can still contribute to the SMA solution of the solar neutrino problem for small 
values of Xe- In addition, another conversion mechanism becomes possible simultaneously if 
Xe ~ Xfi cLS would follow for non-hierarchical Yukawa couplings. The z/g — z/^ mixing angle 
is large in this case, and Ami2 in (^l]) can be appropriate for the vacuum solution when 
xi ~ 10~^ (which also corresponds to the SMA solution). Thus the solar neutrino flux gets 
altered inside the Sun through the resonance with KK states and outside through vacuum 
oscillations with z/^. The simultaneous presence of these two solutions helps in getting better 
agreement with the solar data since it is known that the vacuum oscillations of Ue to sterile 
state gives a poor flt to the rates observed in different solar neutrino experiments. But this 



in conjunction with the SMA MSW resonance with the KK excitations is argued [25| to 



describe various features of the solar neutrino data. However, note that since both the mass 
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scales generated here correspond to Am^, it is not possible to incorporate the solutions to 
both the atmospheric and LSND anomalies through the introduction of a single additional 
neutrino z/^. 

(ii) If n is chosen around the LSND scale, the KK resonance is ineffective as discussed in 
Sec. |T1|. The solar neutrino problem can still find an explanation. This is because the Am^2 
controlling the solar oscillations can now be in the MSW range if /i^ is towards the lower end 
of the LSND region and mixing parameter Xi, ~ 10~^. Since the mixing angle 6efi is large 
(1^), this can provide the large angle MSW solution, which is preferred by the solar data. 
The higher dimensional physics is not directly involved in solar neutrino problem in this 
case but its role is to generate mixing and mass splitting among the active neutrinos. The 
oscillations relevant for the LSND solution could occur indirectly through coupling with u^. 
Incorporating the atmospheric neutrino problem in this context would require introducing 



the third neutrino Ur- This extension is discussed in Sec. [IV C . 

If one or more of the masses fii, fj^2 is zero, the mixing angle generated (^OD is very small. 
Since we need at least one large mixing angle (for solving the atmospheric neutrino prob- 
lem), we shall consider only those scenarios in which at least two neutrinos have a nonzero 
degenerate mass to begin with. If we do not introduce any scale other than the common 
mass fi, we are led to two different possibilities corresponding to (a) two degenerate and one 
massless neutrino (see Sec. [IV B|) , and (b) all three degenerate neutrinos (see Sec. [IV Q ). The 



complete neutrino mass spectrum in these models is determined in terms of five parameters: 
the common mass fi of neutrinos, three mixing parameters C,a and the compactification ra- 
dius R. While their magnitudes are arbitrary, we will concentrate on the consequences that 
follow when they are assumed to be around the following "natural" values: 

^ (2 - 3) • 10^3 eV ; ^« = V2m^R ^ (0.03 - 0.4) ; ^ (g.l - 1) eV^ . (43) 

The value of R is near the observational limit and it also allows the possibility of an MSW 



resonance between u^, and the tower of the KK states [T^. Given this value and assuming 
the fundamental scale Mf to lie in 1 — lOTeV range, one obtains the quoted values for C,a- 
Value of .^e in this range leads to a SMA solution to the solar neutrino problem. Natural 
value of rua ~ lO""^ — 10~^ eV cannot help in generating the LSND or atmospheric mass 
scale. The value of n therefore needs to be chosen near the LSND scale. Given this, the 
atmospheric scale follows naturally as we will see. The solar scale is explained in terms of 
the value of R chosen in the above equation. 



B. Three unmixed generations with Majorana masses {0, ;U,/i} 

We consider the case corresponding to three unmixed neutrinos having the masses 
{0, yU, fi} in the brane. The 4x4 mass matrix involving z/g, u^, Vri is given by 



/ me \ 

yU 

fi 

\me rrir J 



(44) 
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Starting with this tjil and including the corrections due to seesaw approximation as in 
eg ■ (|39D , we get the following effective mass matrix: 



rriL = 



V ^Xe ^Xm ^Xr ) 



(45) 



where r = (y'2E(l/A;2) - EXa/2). 

We can study consequences of the above equation approximately by retaining quadratic 
terms in parameters Xa^ t- This approximation is seen to be quite good for the choice of 
parameters as in eq.(^). The diagonalizing matrix is given in this approximation by 

k = R^m RM) RisW Ruioo) , (46) 

where RijS are 4x4 rotation matrices in the i — j plane, and the angles are given by 

= tan-i ^ , = -TXi. , = ,u; = n/4 (47) 

Xr ^ 



where x± = y xfi + Xt ■ The mass eigenvalues are 

Ai = -fiTXe , Aa = /i , As = /i(l - x±) , A4 = /irxe • (48) 

Note that two of the eigenvalues are degenerate. They will be split by terms higher order 
in x^. By keeping the higher order terms in the diagonalization, one can show that the 1-4 
splitting is of 0{fi^6^) where 6 denotes typical magnitude of Xm ~ Xr ~ Xe ~ t- We thus 
have the following (mass)^ differences: 

Am24 ^ Am34 ^ Aml^ ^ ^^li = A*^ , 

Am?4 ^ C((5V^) • (49) 

The Am23 which plays the role of A^^^ describes the splitting between (almost) degen- 
erate z/^ — Uj. pair. The z/e — pair is almost maximally mixed and is separated from the 
— pail' by the LSND scale. This pattern thus reproduces the 2+2 scheme of neutrino 
mixing. The phenomenology of the solar neutrino however differs considerably from the 2+2 
model. Now the zero mode of vb as well its KK excitations contribute to the solution of 
the solar neutrino problem: the former through non-zero Am^4 and the latter ones through 
MSW resonances with Vf.. This is quantitatively displayed in Fig. 0. This figure shows the 
variation of Aml^ obtained by diagonalizing eq.(^5D with Xe for different values of fi. The 
Xn,T are determined by identifying Am^s and 6 (see eq. ^) with the atmospheric neutrino 
mass scale and mixing angle. It is seen that the favourable value of Xe ~ 0.02 for obtaining 
the SMA solution also corresponds to Am^4 ~ 10~^° — 10~^^eV^. Since the Ue — mixing 
is almost maximal {uj = vr/4), the conversions <-> i^s can take place through vacuum oscil- 
lations. Thus in this model one has the combined effect of the vacuum and SMA resonance 



with the KK states as already discussed in Sec. |IV A 
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The CHOOZ constraint is easily satisfied: The survival probability contains an oscil- 
latory term with the amplitude 

sin^ 2^e ~ ^^l^^l-i (50) 

which vanishes within our approximation. The LSND scale also gives an averaged contribu- 
tion 

which is of (9(10~^) and hence negligible. 

The amplitude of the LSND probability Pipfj, — * is given by^ 

The LSND probability is significantly constrained here by the solutions of the solar and 
atmospheric neutrino problems. The maximum allowed value of ~ X? for which (^Of ) 
reproduces the observed atmospheric scale is approximately given by ~ 0.02. The Xe 
required to be around 3 ■ 10"'* to obtain the SMA solution to the solar neutrino problem. 
As a result, Pe/x < 10~^, which falls short of the LSND observation. 

One can increase the LSND probability in the model by allowing larger value for Xe and 
thus by sacrificing the SMA solution. This is seen from Fig. ^ which shows the effective 
LSND mixing angle following from ( |5TD as a function of Xe- We have chosen Am^j^ and 
Qatm at extreme values in the allowed range so as to maximize X\i. and hence the value of the 
LSND probability in (|^). It is seen from the figure that the observed probability cannot 
be reproduced by the model for Xe < 1 • The trend suggests that one may be able to obtain 
LSND probability if Xe is chosen large. While the SMA solution is no longer there, there is 
an alternative mechanism to solve the solar neutrino problem here. Fig. |^ shows that Am^4 
increases with Xe and eventually for Xe ~ 0.5 — 0.9 one obtains lS.m\^ in the MSW range. 
Thus instead of the higher excitations, the zero mode can cause the MSW transition in this 
case. The perturbative formalism followed here is no longer valid for Xe > 1 and it remains 
to be seen if all neutrino anomalies can be simultaneously understood in this case. 



C. Three unmixed generations with Majorana masses {/i, 

Starting with the mass pattern /z, /x, [i for the active neutrinos, the 4x4 mass matrix 
involving i>e,i'^,i'T, is given by 



rriL 








rrie 





fl 


rrif. 











(52) 



\me. Vflr J 



^ The contributions due to the KK states have been neglected in this approximation. Since from 
(H), the coupling of Ve and to the KK states is ^ • (E ^) ^ {Ue2U^2 + UesU^s) and 
decreases rapidly (oc 1/n^) with increasing n, this approximation is valid up to a factor of 0(1). 
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Starting with this ttil and including the corrections due to seesaw approximation as in 
eg ■ (|39D , the effective mass matrix is obtained as 



fi-xl 



'XeXfJ. XeXr TXt 



\ 



XeXfJ. 1 Xfj, XfiXr TX/i 
XeXr XfiXr 1 Xt ^Xt 

V -rXe rxt, tXt / 



(53) 



The above matrix can be diagonahzed exactly. The eigenvalues are given by 

A, = /X , = /X , A3,4 = ^ ((1 - X') ± [(1 - xr + 4r^xT/') , 
where x"^ = xl + xfi + Xt- This leads to the following (mass)^ differences: 



(54) 



Amf2 = 



Amis 



2^,2 



2 

23 ~ ^1^ A. 5 
2 _ A™2 . . A_2 _ ,.2 



The diagonalizing matrix is given up to terms quadratic in {xa, t} by 

K = R23{e) R3^{<P) RM) Ruioo) , 



where the angles are given by 



6 = tan 



Xt 



—TX± , 'ip = tan 



Xe_ 

x± 



TXeXl. 

id 



(55) 



(56) 



(57) 



This model contains two exactly degenerate states. As a result, one has only two independent 
Am^s and it is not possible to account for all the anomalies. Moreover, the mixing pattern 
in ( |56D is such that even if one is willing to give up LSND, the solar and atmospheric 
neutrino anomalies cannot be simultaneously explained. In order to do this, one would need 
to identify the larger Aw? ^ fi with the atmospheric neutrino scale. But in that case, 
the atmospheric neutrino mixing angle following from (^) turns out be too small [(9(s^)]. 
In spite of vanishing Am^2; the solar anomaly can be accounted through resonance with 
the KK states. As shown in section (3), this becomes feasible only if the electron neutrino 
mass fi < 10~^eV. In this case, the Am^s in eq. (p5|) cannot account for the atmospheric 
neutrino anomaly. Thus simultaneous explanation of the solar and atmospheric neutrino is 
not possible and the model does not seem phenomenologically viable. 



V. ZEE MODEL FOR SOLVING ALL ANOMALIES 



The neutrinos in the brane were assumed unmixed and degenerate so far. Neutrino 
oscillations occurred entirely due to the presence of coupling to the bulk neutrino. Since 
that seems to be inadequate for explaining all the neutrino anomalies, we now consider a 
more general possibility in which some of the brane neutrino mixings are present even in 
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the absence of the bulk states. The latter can provide additional structure needed to under- 
stand all neutrino anomalies. This possibility will make some of the simple neutrino mass 
generation mechanisms viable which by themselves cannot solve all the neutrino anomalies. 
This is exemplified by the model due to Zee [^. We confine our discussion to the Zee model 



although the basic formalism in Sec. y can be used for any arbitrary mass structure in the 
brane. 

One obtains the following neutrino mass matrix in the Zee model: 

/ e se\ 



e C0 
\S0 Cg J 



where 



tan 6' = 



f, 



e/1 



(58) 



(59) 



The /a/3 in the above equation are Yukawa couplings of the charged singlet Higgs to the 
leptonic doublet and /i is the overall mass scale. 

The above structure has been used to simultaneously solve the solar and atmospheric 
p7| or the LSND and the atmospheric neutrino anomalies pHj^. The non-hierarchical / 
imply a very small e in eq. ( ^81) , leading to an approximate + — symmetry. This 
corresponds to maximally mixed degenerate pairs. Only — z/g oscillations occur in this 
limit and the LSND result can be explained by choosing 4sg 10"'^ and /i^ ~ 0.1 — 1 eV^. 
Non-zero e can generate spitting of O^i^'^esg) between the degenerate pairs. This would 
correspond to the atmospheric neutrino scale if e ~ 10~^. While small values for e are 
more natural in this model, the required magnitudes of e, sg can be possible with somewhat 
inverted hierarchy f^r ^ /^t ^ /e/^- The solar neutrino problem cannot be accommodated in 
the Zee model with this choice of parameters. This however becomes possible once coupling 
to bulk neutrino is switched on. This coupling also allows generation of the atmospheric 
neutrino scale even when e is zero, i.e. model is Lg + — Lr symmetric. 

The coupling of Zee model to the bulk neutrino leads to the following mass matrix mi 
in the basis Ueyi'fj,, Hr, to zeroth order in the seesaw approximation: 



/ 



rriL = fJ, 



Sg r, 
t 

Sg Cg Tr 



\ 



Te 



(60) 



) 



Including non-leading corrections, we have the following effective mass matrix for the four 
neutrino states: 

5mT= ^ 



X 

2 

XrXe^, + eYe/. 2(ceX^Xr + tXeX^l) 

SgYer + CgXeXf^ + eX/^Xr CgY^^ + SgXeX^l + (^XeXr 



SgYer + CgXeX^i + eX/xXr Xef3 \ 

CgY^r + SgXeX^i + eXeXr X^l^^ 

'^iXrXeSg + XrX^iCg) XrP 

XrP / 



(61) 
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where 



a/3 



Ta = main , P = iXeVe + Xm^M + Xr^r) • 



(62) 



Note that for the natural values m^, ~ 10~^ eV and for fi ~ eV, the terms proportional 
to f3 are much smaller than other elements in the matrix and can be neglected. Then up to 
the second order in other parameters, tjil can be written as: 



rriL 



( e - XeXr/2 Sq - XeXM/2 \ 

e - XeXr/2 -x^Xr 1 - 8^2 - xi/2 

SB - XeX^/S 1 - sll2 - xi/2 -x^Xr 

V oy 



(63) 



This matrix can be diagonalized through 



R 



,se - e 



XeX- 



2 ' "''^ V2 2 ' "^•'^ 4 

where x± = {x^i i Xt) / V^- The rotation matrix (X) may be expanded as 



■) R23{ 



p2 _ „2 
t — bn 



) 



(64) 



2 



71^ 
1 

V2 



XeX- 



„2 

2 2 ' 



1 



-1 



_ XeX + 

V2 ^ 2 

'1 _ £s 

^ 2 

_2 



V 



j_ £ _ J_('_1 _L ^ 

^2 2 ) ^\ ^^2 





He ' 
2 ' 







^) 



(65) 



ly 

The mass eigenvalues are 

Ai = -2/iese , A2 = //(-I + X- + - eV2) , A3 = //(I - X+ + + eV2) , A4 - 0. (66) 
Thus, the three mass squared differences are 



Am24 ^ Am34 ^ Am2i ^ Amg^ fsi /^^ , 
Arr4 ^ 4/i2(x^XT + ese) , 
Ami4 pa 4yLi^e^s^ . 



(67) 



The mass pattern is similar to the conventional 2+2 scheme. This allows the solution to all 
anomalies simultaneously: The amplitude of the LSND oscillations is given by 



(68) 



Thus, the choice of sq made at the beginning of the section is a suitable one. All that is 
needed is = Arnj^g^jj. The large mixing required for the atmospheric mixing is naturally 
obtained: here 



(69) 
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so that the mixing is nearly maximal, and the mass squared difference is 



4/i'(XMXr + ese) = ^ml,^ ^ 3 x 10'=^ eV^ . (70) 

Note that ( pE| ) and ( pPD do not constrain the value of e in any way, but (^) restricts the 
maximum value that e can take: if the major contribution to Am^^^ is from the e term, we 
have e ~ Am^^^/(4/i^S6i) ~ 10~^ — 10~^. The z/g survival probability contains an oscillatory 
term with the amplitude 

sin^ 2^, ^ ^\Jl,\Jl ^ {si - e'f < lO'^ . (71) 

The LSND scale also adds an average term to this probability with the amplitude 2(e^ + s^). 
Both these amplitudes are within the observed bounds in the CHOOZ experiment. Notice 
that this experiment puts a similar upper bound on the allowed value of e as the atmospheric 

For e ~ 0.1 and se ~ 10^^, the smallest mass difference Amf^ ~ 10"^ - lO^^eVl While 
this is the right value for the MSW effect, the resonance cannot occur with the massless 
mode due to the fact that the electron neutrino is heavier. The v^. can however resonate 
with the KK excitations with mass ^ > Ai. Thus the solar neutrino problem can be solved 
as in il^, with quantitative details differing due to the presence of non-zero mass Ai. For 



e ~ 0.1 and 6 ~ 0.01, this mass is in the MSW range and this can significantly alter the 



survival probability in the manner discussed in Sec. [ITT . 

Let us now consider the limiting case of the Zee model obtained when e ^ 0. As already 
mentioned, the mass matrix is invariant under + — Lj. in this case and implies a 
degenerate — ^'/i pair- Such a matrix cannot lead to the atmospheric neutrino scale. 
This scale is generated through the couplings to bulk neutrinos which violate Le + L^ — Lt 
symmetry. The expressions for the mixing and masses can be recovered from the earlier case 
by putting e = 0. This limit does not affect the solution to neutrino anomalies. One can 
simultaneously solve all anomalies for the similar values of parameters i?, sq, h, Xa as in the 
case with non-zero e. The electron neutrino mass now is much smaller than the MSW scale. 



As a result, the solution to solar neutrino problem occurs exactly as in [|T7| with very little 
perturbation from Ai. 



VI. SUMMARY 

We have analyzed the neutrino mass spectra in models based on higher dimensional 
theories with an extra dimension of mm size. We have restricted our discussion to the 
minimal scenario where a single fermion propagating in the bulk couples weakly to the 
flavour neutrinos in the brane. This coupling can generate and / or modify the masses and 
mixings in flavour neutrinos. In addition, the singlet fermion in the bulk provides a massless 
sterile neutrino and a KK tower of several light sterile states which can contribute to neutrino 
oscillations. We have developed a formalism for calculating the net neutrino mass spectrum, 
taking into account the effect of the brane-bulk coupling. Simple approximations allow one 
to discuss various features of the neutrino mass spectrum analytically. We examined the 
possibility of solving all the neutrino anomalies — atmospheric, solar and LSND — through 
this in a natural manner. 
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We have calculated the masses and mixings for several cases that can potentially solve 
all the neutrino anomalies. In the absence of any masses for the brane neutrinos, it is 
not possible to generate the required masses and large mixings naturally just through the 
coupling with the bulk fermion. If two or more neutrinos are massive and degenerate, 
however, large mixings and hierarchical mass splittings are possible. We considered two 
cases in detail, the one in which all three flavour neutrinos are degenerate, and the one in 
which two of them are degenerate and the third massless. It turns out that if all the three 
neutrinos are degenerate, their coupling to the bulk fermion cannot lead to a simultaneous 
understanding of even the solar and atmospheric neutrino anomalies, let alone LSND in 
addition. This leaves the other option as the only viable alternative. 

When two of the three neutrinos are degenerate with mass /i and the neutrinos have 
no a priori mixings among themselves, their couplings to the bulk fermion automatically 
lead to two pairs of almost degenerate neutrinos separated by the scale characterized by /i. 
The splittings within the pair are hierarchical, and may account for Am^ and Am^^^. The 
hierarchy is completely controlled by the higher dimensional physics and one typically finds 
Amg/Am^j^ ~ 5^ where 6 refers to the typical mixing between the bulk and brane neutrinos. 
This hierarchy generates the scale corresponding to long wavelength solar oscillations for 5"^ ~ 
10~^. Identification of fi with the LSND scale turns out to generate the atmospheric scale 
for the same value of 6. Thus the features needed to account for the solar and atmospheric 
neutrino anomalies follow automatically. But the z/e ^ z/^ oscillation probability is much 
smaller than that observed at LSND when the mixing 6 is small. 

Understanding of all neutrino anomalies would then require some preexisting mass struc- 
ture among the flavour neutrinos. We analyzed the specific example of the neutrino mass 
matrix predicted by the Zee model, which produces two massive degenerate neutrinos and a 
massless one. The coupling with the bulk fermion leads to the 2+2 mass spectrum, and its 
solutions to solar and atmospheric neutrino anomalies can coexist with large Ue — z/^ oscilla- 
tion probability seen at LSND. Even in the e ^ limit of the Zee model (which corresponds 
to the Le + L^ — Lr symmetry), all the three neutrino anomalies can be naturally accounted 
for. The LSND scale /i, the LSND mixing angle 6 and the large atmospheric mixing angle 
are already present in the structure of the Zee (or + — Lj. symmetric) mass matrix. 
The bulk fermion provides a massless sterile neutrino and a KK tower of sterile states that 
can participate in the neutrino oscillations. The size of extra dimensions R ~ mm creates 
masses of the lighter KK states in the range m\ ~ 10~^ eV^, appropriate for the SMA 
MSW solution of the solar neutrinos. The brane-bulk coupling also generates the splitting 
^^atm ~ and the solar mixing angle sin^Q ~ ^e/n. The addition of the bulk fermion 

and its coupling simultaneously provides Am^^^, Am| and 6q in the right range, thus com- 
pleting the picture. The Zee model embedded in extra dimensions in the minimal way can 
then account for all the neutrino anomalies. 

The sterile neutrinos participate in the solar neutrino oscillations. We have calculated 
how the survival probabilities in get modified due to the nonzero masses of z/g that is 
obtained in the preferred scenarios. Another interesting feature of some of the scenarios 
considered here is the simultaneous occurrence of two different mechanisms for the solution 
to the solar neutrino problem. The oscillations of z/g to massless mode of the bulk state 
corresponds to the long wavelength solution and oscillations to higher modes are appropriate 
for the SMA MSW conversion. These lead to differences compared to the phenomenology of 
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the conventional "2+2" schemes [|TT| and allows one to explain ||25| various features of the 
solar neutrino data. 

We had assumed throughout that size of the extra dimension and the fundamental scale 
are near their observable limits. It is seen from the present considerations that this ob- 
servability does not conflict with the observed features of neutrino masses and mixings. 
While the presence of large extra dimensions cannot exclusively account for all the neutrino 
anomalies they can provide an important ingredient for generating the observed features of 
the neutrino spectrum. 
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FIGURES 



surv 1.0 




E (MeV) 

FIG. 1. Survival probability of from the sun as a function of energy. The (solid, long dashed, 
short dashed, long-short dashed, long-short-short dashed) curves are for A = (0, l,2,5,103)-10-3eV. 
The purpose of this plot is only to illustrate the dependence of Pgurv on A. The details of neutrino 
production and transport inside the sun and the energy resolution have not been taken into account. 
These will result in a smearing of the energy dependence shown in the figure. 
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FIG. 2. Am|i in eV^ (the lower pair of curves) and the predicted LSND probability (eq. |5l[) 



^"'^ (the middle pair of curves) 



,2 



plotted against Xe- The solid (dotted) curves correspond to 
= 1.0 (0.2) eV^. The uppermost horizontal lines are experimental LSND probabilities .02 and 
.003 corresponding to jj? = 1.0 (solid) and 0.2 eV^ (dotted) respectively. 
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